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Abstract 

We apply the mathematical technique of factorization of differen- 
tial operators to two different problems. First we review our results 
related to the supersymmetry of the Montroll kinks moving onto the 
microtubule walls as well as mentioning the sine-Gordon model for the 
microtubule nonlinear excitations. Second, we find analytic expressions 
for a class of one-parameter solutions of a sort of diffusion equation of 
Bessel type that is obtained by supersymmetry from the homogeneous 
form of a simple damped wave equations derived in the works of P.A. 
Robinson and collaborators for the corticothalamic system. We also 
present a possible interpretation of the diffusion equation in the brain 
context. 



1 Nonlinear biological excitations 

The possibility of soliton excitations in biological structures has been first 
pointed out by Englander et al Q.] in 1980 who speculated that the so-called 
'open states' units made of approximately ten adjacent open pairs in long 
polynucleotide double helices could be thermally induced solitons of the 
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double helix due to a coherence of the twist deformation energy. Since then 
a substantial amount of literature has been accumulating on the biological 
significance of DNA nonlinear excitations (for a recent paper, see [2]). On 
the other hand, the idea of nonlinear excitations has emerged in 1993 in 
the context of the microtubules (MTs) [3J, the dimeric tubular polymers 
that contribute the main part of the eukaryotic cytoskeleton. In the case 
of neurons, MTs are critical for the growth and maintenance of axons. It 
is known that axonal MTs are spatially organized but are not under the 
influence of a MT-organizing center as in other cells. We also remind that 
in 1995 Das and Schwarz have used a two-dimensional smectic liquid crystal 
model to show the possibility of electrical solitary wave propagation in cell 
membranes [3]. Nevertheless, there is no clear experimental evidence at the 
moment of any of these biological solitons and kinks. 



2 Supersymmetric MT Kinks 

Based on well-established results of Collins, Blumen, Currie and Ross [5] re- 
garding the dynamics of domain walls in ferrodistortive materials, Tuszyfiski 
and collaborators [3j [6] considered MTs to be ferrodistortive and studied 
kinks of the Montroll type [Zj as excitations responsible for the energy trans- 
fer within this highly interesting biological context. 

The Euler-Lagrange dimensionless equation of motion of ferrodistortive 
domain walls as derived in [5] from a Ginzburg-Landau free energy with 
driven field and dissipation included is of the travelling reaction-diffusion 
type 

iff" + pip' - V 3 + V> + ° = , (!) 

where the primes are derivatives with respect to a travelling coordinate 
£ = x — vt, p is a friction coefficient and a is related to the driven field [5]. 

There may be ferrodistortive domain walls that can be identified with 
the Montroll kink solution of Eq. (fTJ) 

where (3 = (pt2 — ct\)/\/2 and the parameters ot\ and a 2 are two nonequal 
solutions of the cubic equation 

(ip - ai)(ip - a 2 ){ip - a 3 ) = ip 3 - tp - o . (3) 
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Rosu has noted that MontrolFs kink can be written as a typical tanh 
kink [8] 

M(£)=7-tanh(^) , (4) 

where 7 = a\ + a<i = 1 + ■ The latter relationship allows one to use 
a simple construction method of exactly soluble double-well potentials in 
the Schrodinger equation proposed by Caticha [9]. The scheme is a non- 
standard application of Witten's supersymmetric quantum mechanics |10j 
having as the essential assumption the idea of considering the M kink as 
the switching function between the two lowest eigenstates of the Schrodinger 
equation with a double-well potential. Thus 

01 = M<p , (5) 

where </>o,i are solutions of (f)' Q1 + [eo,i — ^(£)]0o,i(£) = 0, and is the 
double-well potential to be found. 
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Figure 1: Single electron within the traveling double- well potential u(£) as a qubit. 
The electron can switch from one wall to another by tunneling and the relation 
between the wavefunctions in the two wells is given by Eq. 

Substituting Eq. ([5]) into the Schrodinger equation for the subscript 1 
and substracting the same equation multiplied by the switching function for 
the subscript 0, one obtains 

0o + Rm^o = , (6) 
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where Rm is given by 



Rri 



M + eM 
2M 1 



(7) 



and e = ei — eo is the lowest energy splitting in the double-well Schrodinger 
equation. In addition, notice that Eq. ([6]) is the basic equation introducing 
the superpotential R in Witten's supersymmetric quantum mechanics, i.e., 
the Riccati solution. For Montroll's kink the corresponding Riccati solution 
reads 



Rm(0 
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tanh ( -£ I + 



2(3 



sinhQ3£) + 2 7 cosh 2 ( -£ 



(8) 



and the ground-state Schrodinger function is found by means of Eq. ([6]) 



0o,m(C) 



0o(O) cosh ( ^£ ) exp ( ) exp ( — — ^ cosh(/3 £) 



2/5 2 



- 7 /3£-7sinh(/3£) 



0) 



while is obtained by switching the ground-state wave function by means 
of M. This ground-state wave function is of supersymmetric type 



</>0,m(£) = </>o,M(0)exp 



Ru{y)dy 



(10) 



where </>o,m(0) is a normalization constant. 

The Montroll double well potential is determined up to the additive 
constant eo by the 'bosonic' Riccati equation 

«Af(c) = - #M + eo = -r + + ^ + £ o 



+- 



(4 7 2 e + 2( 7 2 + l)ecosh(/3£) - 8/3 2 ) cosh(/3£) 

-4 7 (e + ecosh(/3£) - 2/3 2 ) sinh(/3f ) 
If, as suggested by Caticha, one chooses the ground state energy to be 



fill 
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e e 
2 + 4^2 



(i-y), 



(12) 



then «m(0 turns into a travelling, asymmetric Morse double- well potential 
of depths depending on the Montroll parameters and 7 and the splitting 
e 



U, 



L,R 



0,m 



P 2 



1 ± 



2e 7 



(13) 
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where the subscript m stands for Morse and the superscripts L and R for 
left and right well, respectively. The difference in depth, the bias, is A m = 



Uq = 2e7, while the location of the potential minima on the traveling 



axis is at 



-_L,R 



(14) 



that shows that 7 / ±1. 



An extension of the previous results is possible if one notices that Rm 
in Eq. ([8]) is only the particular solution of Eq. (jlll) . The general solution 
is a one-parameter function of the form 



R M {t; A) = R M (0 + -| [ln(I M (0 + A) 



and the corresponding one-parameter Montroll potential is given by 

cP_ r. 



um(0~^ ln(/M(6 + A) 



(15) 



(16) 



In these formulas, Im(0 = 0o m(£)^£ an d A is an integration constant 
that is used as a deforming parameter of the potential and is related to the 
irregular zero mode. The one-parameter Darboux-deformed ground state 
wave function can be shown to be 



(£;A) = J\(\ + l 



<Po,M 



Im(C) + A 



(17) 



where \J A(A + 1) is the normalization factor implying that A ^ [0,-1]. 
Moreover, the one-parameter potentials and wave functions display singular- 
ities at A s = — Im(£,s)- For large values of ±A the singularity moves towards 
=f 00 and the potential and ground state wave function recover the shapes of 
the non-parametric potential and wave function. The one-parameter Morse 
case corresponds formally to the change of subscript M — > m in Eqs. (|15p 
and (I16p . For the single well Morse potential the one-parameter procedure 
has been studied by Filho [12] and Bentaiba et al [15] . 

The one-parameter extension leads to singularities in the double- well po- 
tential and the corresponding wave functions. If the parameter A is positive 
the singularity is to be found on the negative f; axis, while for negative A it is 
on the positive side. Potentials and wave functions with singularities are not 
so strange as it seems |14j and could be quite relevant even in nanotechnol- 
ogy where quantum singular interactions of the contact type are appropriate 
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for describing nanoscale quantum devices. We interpret the singularity as 
representing the effect of an impurity moving along the MT in one direc- 
tion or the other depending on the sign of the parameter A. The impurity 
may represent a protein attached to the MT or a structural discontinuity in 
the arrangement of the tubulin molecules. This interpretation of impurities 
has been given by Trpisova and Tuszyhski in non-supersymmetric models of 
nonlinear MT excitations [15] . 

3 The sine-Gordon MT solitons 

Almost simultaneously with Sataric, Tuszynski and Zakula, there was an- 
other group, Chou, Zhang and Maggiora [16], who published a paper on 
the possibility of kinklike excitations of sine-Gordon type in MTs but in a 
biological journal. Even more, they assumed that the kink is excited by the 
energy released in the hydrolysis of GTP — ► GDP in microtubular solutions. 
As the kink moves forward, the individual tubulin molecules involved in the 
kink undergo motion that can be likened to the dislocation of atoms within 
the crystal lattice. 

They performed an energy estimation showing that a kink in the system 
possesses about 0.36 - 0.44 eV, which is quite close to the 0.49 eV of energy 
released from the hydrolysis of GTP. 

Moreover, they assumed that the interaction energy U(r) between two 
neighboring tubulin molecules along a protofilament is harmonic: 



where k = — -A °^ and r = x% — In addition to this kind of nearest 

neighbor interaction, a tubulin molecule is also subjected to interactions 
with the remaining tubulin molecules of the MT, i.e., those in the same 
protofilament but not nearest neighbor to it. 

Chou et al cite pages 425-427 in the book of R.K. Dodd et al (Soli- 
tons and Nonlinear Wave Equations, Academic Press 1982) for the claiming 
that this interaction for the ith tubulin molecule of a protofilament can be 
approximated by the following periodic effective potential 



where XJq is the half-height of the potential energy barrier and £j is the dis- 
placement of the ith tubulin molecule from the equilibrium position within 
a particular protofilament. 



U(r) « \k{r-a f 



(18) 




(19) 



6 



Introducing the new variable fa = §^£j the following sine-Gordon equa- 
tion is obtained 



d 2 (j) , 2 d 2 <P /2tt\ 2 . , 
m 1^2 = ka o*Z2 ~ — ^oW (20) 



dt 2 "dx 2 \a J 
that can be reduced to the standard form of the sine-Gordon equation 

d 2 <p 1 d 2 4> 1 



dx 2 c 2 dt 2 I 2 



sin</> (21) 



if one sets c 2 = ^ and I 2 = 4 ^ 2 ^ . Now, it is well known that the 
sine-Gordon equation has the famous inverse tangent kink solution 

= tan" 1 (exp[±j(x-vt)]j , (22) 

where 7 = — , 1 is an acoustic Lorentz factor and w = ? is the kink width. 

Most interestingly, the momentum of a tubulin dimer is strongly local- 
ized: 

p = ——. — - = — sech — -(x — vt) . (23) 

at 7r I [ I 

This momentum function possesses a very high and narrow peak at the 
center of the kink width implying that the corresponding tubulin molecule 
will have maximum momentum when it is at the top of the periodic potential. 
According to Chou et al this remarkable feature occurs only in nonlinear 
wave mechanics. 

Interestingly, for purposes of illustration, these authors have assumed 
the width of a kink w ~ 3ao- Therefore, with the kink moving forward, the 
affected region always involves three tubulin molecules. For a general case, 
however, the width w of a kink can be calculated from 




if the force constant k between two neighboring tubulin molecules along a 
protofilament, the distance ao of their centers, and the energy barrier 2Uq 
of the periodic, effective potential are known. Then the number of tubulin 
molecules involved in a kink is given by 

- = (2^~ 1 Jka 2 /U Q . (25) 
ao v 
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It is further known that the tubulin molecules in a MT are held by nonco- 
valent bonds, therefore the interaction among them might involve hydrogen 
bonds, van der Waals contact, salt bridges, and hydrophobic interactions. 

It was found by Israelachvili and Pashley [T7] that the hydrophobic force 
law over the distance range 0-10 nm at 21°C is well described by 

EE. = Ce~ D / Do N/m , (26) 
R 

where D is the distance between tubulin molecules, Dq is a decay length, and 
R = R1+R2 * s a harmonic mean radius for two hydrophobic solute molecules, 
all in nm. R is 4 nm in the case of tubulin. 

3.1 More on the hydrolysis and solitary waves in MTs 

Inside the cell, the MTs exist in an unstable dynamic state characterized 
by a continuous addition and dissociation of the molecules of tubulin. The 
polypeptides a and (5 tubulin each bind one molecule of guanine nucleotide 
with high affinity. The nucleotide binding site on a tubulin binds GTP 
nonexchangeably and is referred to as the N site. The binding site on (3 
tubulin exchanges rapidly with free nucleotide in the tubulin heterodimer 
and is referred to as the E site. 

Thus, the addition of each tubulin is accompanied by the hydrolysis of 
GTP 5' bound to the (3 monomer. In this reaction an amount of energy of 
6.25 x 10~ 21 J is freed that can travel along MTs as a kinklike solitary wave. 

The exchangeable GTP hydrolyses very soon after the tubulin binds to 
the MT. At pH = 7 this reaction takes place according to the formula: 

GTP 4 - + H 2 -> GDP 3 - + HPOj- + H + + A H E . (27) 

The last mathematical formulation of the manner in which the energy 
AhE is turned into a kink excitation claims that the hydrolysis causes a 
dynamical transition in the structure of tubulin |18| . 

4 Quantum information in the MT walls 

Biological information processing, storage, and transduction occurring by 
computer-like transfer and resonance among the dimer units of MTs have 
been first suggested by Hamerrof and Watt [19] and enjoys much speculative 
activity [20] . 

For the case of sine-Gordon solitons, the information transport has been 
investigated by Abdalla et al [2~T] . 
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Recently Shi and collaborators [22] worked out a processing scheme of 
quantum information along the MT walls by using previous hints of Lloyd 
for two- level pseudospin systems |23| . The MT wall is treated as a chain of 
three types of two pseudospin-state dimers. A set of appropriate resonant 
frequencies has been given. They conclude that specific frequencies of laser 
pulse excitations can be applied in order to generate quantum information 
processing. 

Lloyd's scheme uses the driving of a quantum computer by means of a 
sequence of laser pulses. He assumes a 1-dimensional arrangement of atoms 
of two types (A and B) that could be each of them in one of two states and 
are affected only by nearest neighbors. Then, information processing could 
be performed by laser pulses of specific frequencies 0JK af3 , that change the 
state of the atom of the K kind (A or B type) if in a pair of atoms AB, A 
is in a state and B is in state (3. 

5 Supersymmetry at the Brain Scale 

Neuronal activity is the result of the propagation of impulses generated at 
the neuron cell body and transmitted along axons to other neurons. Re- 
cently, Robinson and collaborators [25] obtained simple damped wave equa- 
tions for the axonal pulse fields propagating at speed v a between two popu- 
lations, a and b, of neurons in the thalamocortical region of the brain. The 
explicit form of their equation is 

d R Mt) = S[V a (t)] , (28) 

where 

A /id 2 2d 9 9 \ 

" = {^ + ^Jt + 1 - r ' v ) • < 29) 

where v a = v a /r a , r a is the mean range of axons a, and V a = J2b^ab is a 
so-called cell body potential which results from the filtered dendritic tree 
inputs. Robinson has used the experimental parameters in this equation for 
the processing of the experimental data. In the following we concentrate on 
a particular mathematical aspect of this equation and refer the reader to 
the works of Robinson's group for more details concerning this equation. 

5.1 The homogeneous equation 

We treat first the homogeneous case, i.e., 5 = and we discard the subindexes 
as being related to the phenomenology not to the mathematics. Let us em- 
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ploy the change of variable z = ax+by—ct (see, e.g., |26j). which is a traveling 
coordinate in 2+1 dimensions. This is justified because it was noticed by 
Wilson and Cowan [24] that distinct anatomical regions of cerebral cortex 
and of thalamic nuclei are functionally two-dimensional although extending 
to three spatial coordinates is trivial. We have the following rescalings of 
functions: <j) t = -ccp z , <j) a = v 2 (j) zz , <j) xx = a 2 (j) zz , <j) yy = b 2 (p zz . Then, we 
get the ordinary differential equation corresponding to the damped wave 
equation in the following form 

6li ' z(t) = (£2 ~ 2fl lz~ + ^) = ' (30) 

where 



vc 



^ c 2 — u 2 r 2 (a 2 + b 2 ) ' [c 2 — u 2 r 2 (a 2 + 6 2 )] 2 ^ 

The simple damped oscillator equation ([301) can be easily factorized 

lu=(4--A(4--A<i>=^<i>- ( 32 ) 



dz J \dz 

The case c 2 < v 2 r 2 (a 2 + b 2 ) implies fj, < and the general solution of (|30p 
can be written 

(j)(z) = e^ z {Ae az + Be~ az ) . (33) 

The opposite case c 2 > z^ 2 r 2 (a 2 + b 2 ) will lead to only a change of sign in 
front of fi in all formulas henceforth, whereas the case c 2 = z^ 2 r 2 (a 2 + b 2 ) 
will be considered as nonphysical. The non-uniqueness of the factorization 
of second-order differential operators has been exploited in a previous paper 
[27] on the example of the Newton classical damped oscillator, i.e., 

^(J^i^)'"- ' <34) 

which is similar to the equation (|30p . unless the coefficient 2(3 is the friction 
constant per unit mass, loq is the natural frequency of the oscillator, and 
the independent variable is just time not the traveling variable. Proceed- 
ing along the lines of [27] , one can search for the most general isospectral 
factorization 

(D z + f(z))(D z +g(z))<i ) = a 2 ( j ) . (35) 
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After simple algebraic manipulations one finds the conditions / + g = — 2fi 
and dg/dz + fg = /x 2 having as general solution fx = x ^ +1 — fi, whereas 
/o = — (jl is only a particular solution. Using the general solution f\ we get 

i +A i_ A = (d z + — ^ - fi \ (d z - - m) 4> = o?<f> . (36) 



This equation does not provide anything new since it is just equation (|3ip. 
However, a different operator, which is a supersymmetric partner of (j36[) 
is obtained by applying the factorizing A-dependent operators in reversed 
order 

A^A + J s (d.-^-^(d. + j JL--^ a > } . (37) 
The latter equation can be written as follows 



or 



where 



£-n4 +w > w )*=o, m 



U \z) = „ 2 - a? - -±—1 (40) 

is a sort of parametric angular frequency with respect to the traveling coor- 
dinate. 

This new second-order linear damping equation contains the additional 
last term with respect to its initial partner, which may be thought of as 
the Darboux transform part of the frequency [28]. Z\ = 1/A occurs as a 
new traveling scale in the damped wave problem and acts as a modulation 
scale. If this traveling scale is infinite, the ordinary damped wave problem 
is recovered. The <f> modes can be obtained from the 4> modes by operatorial 
means |27| . 

Eliminating the first derivative term in the parametric damped oscillator 
equation (|39|) one can get the following Bessel equation 

d 2 u ( v? — \ „\ , 
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where x = z + 1/A, n 2 = 5/4, and = ia. Using the latter equation, the 
general solution of equation (|39[) can be written in terms of the modified 
Bessel functions 

4> = (z + l/Xf/^CJ^aiz + 1/A)) + C 2 I_^ /2 (a{z + 1/A))]e^ . (42) 

What could be a right interpretation of the supersymmetric partner equa- 
tion ([37]) ? Since the solutions are modified Bessel functions, we consider 
this equation as a diffusion equation with a diffusion coefficient depending 
on the traveling coordinate. Noticing that the velocity in the traveling vari- 
able of this diffusion is the same as the velocity of the neuronal pulses we 
identify it with the diffusion of various molecules, mostly hormones, in the 
extracellular space (ECS) of the brain, which is known to be necessary for 
chemical signaling and for neurons and glia to access nutrients and thera- 
peutics occupying as much as 20 % of total brain volume in vivo [29J. 

5.2 The nonhomogeneous equation 

The source term S in Robinson's equation (|28p is a sigmoidal firing function, 
which despite corresponding to a realistic case led him to work out extensive 
numerical analyses. Analytic results have been obtained recently by Troy 
and Shusterman [30J by using a source term comprising a combination of 
discontinuous exponential coupling rate functions and Heaviside firing rate 
functions. In addition, Brackley and Turner [31J incorporated fluctuating 
firing thresholds about a mean value as a source of noisy behavior |31j . 

The procedure of Troy and Shusterman can be applied for the paramet- 
ric damped oscillator equation as well as for the Bessel diffusion equation 
obtained herein in the realm of Robinson's brain wave equation with the 
difference that the method of variation of parameters should be employed. 
The detailed mathematical analysis is left for a future work. 
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